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Abstract 

We comment on theoretical aspects of the measurement of azimuthal asymmetries 
in semi-inclusive charged particle production, made recently by the ZEUS Collabo- 
ration at HERA. By taking the ratio of the two measured asymmetries, we find good 
agreement between the perturbative QCD prediction and the experimental data. To 
separate the perturbative and nonperturbative contributions to the asymmetries, we 
suggest that the azimuthal asymmetries of the transverse energy flow be measured 
as a function of a variable qr related to the pseudorapidity of the energy flow. 
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In a recent publication [1] the ZEUS Collaboration at DESY-HERA has pre- 
sented data on asymmetries of charged particle (h^) production in the process 

7* 4- 

e+p — ► e + hr 1 + X, with respect to the angle tp defined as the angle between 
the lepton scattering plane and the hadron production plane (of and the 
exchanged virtual photon). The azimuthal asymmetries, (cosy) and (cos2y), 
as functions of the minimal transverse momentum p c of the observed charged 
hadron in the hadron-photon center-of-mass (hCM) frame, are defined as 

(cosny)(p c ) = Jd y° d ^ COSUip ^dzdQ'M PTd ^ (1) 

with n = 1,2. In terms of the momenta of the initial proton P M , the final- 
state hadron P%, and the exchanged photon g M , the variables in (1) are Q 2 = 
-q^, x = Q 2 /2(P ■ q), and z — (P • Ph)/(P • q). fd$ denotes the integral 
over x,z,Q 2 ,Pt within the region defined by 0.01 < x < 0.1, 180 GeV 2 < 
Q 2 < 7220 GeV 2 , 0.2 < z < 1, and p T > p c - Nonzero (cos 2<p) comes from 
interference of the helicity +1 and —1 amplitudes of the transverse photon 
polarization; and nonzero (cosy) comes from interference of transverse and 
longitudinal photon polarization. 

More than 20 years ago it was proposed to test QCD by comparing measured 
azimuthal asymmetries to the perturbative predictions [2]. However, it was also 
realized that nonperturbative contributions and higher-twist effects may affect 
the comparison [3-6]. For example, intrinsic kr might be used to parametrize 
the nonperturbative effects [3], and indeed ZEUS did apply this idea to their 
analysis of the data [1]. The relative importance of the nonperturbative effects 
is expected to decrease as pr increases. Thus, the azimuthal asymmetries in 
semi-inclusive deep-inelastic scattering (SIDIS) events with large pr should 
be accurately described by perturbative dynamics. From the comparison to 
the perturbative QCD calculation at the leading order in a s [7,8], the ZEUS 
Collaboration concluded that the data on the azimuthal asymmetries at large 
values of p c , although not well described by the QCD predictions, do provide 
clear evidence for a perturbative QCD contribution to the azimuthal asym- 
metries. 



In this paper, we will take a new look at the ZEUS data, motivated by a 
QCD resummation formalism [9-12] that takes into account the effects of 
multiple soft parton emission. First, we argue that the analysis of (cosy?) 
and (cos2yj) based on fixed-order QCD is unsatisfactory because it ignores 
large logarithmic corrections due to soft parton emission. We also show that 
perturbative and nonperturbative contributions are mixed in the transverse 
momentum distributions. Then we make two suggestions for improvement of 
the analysis of the ZEUS data. We show that perturbative and nonperturbative 
contributions can be separated more clearly in asymmetries depending on a 
variable qr related to the pseudorapidity of the final hadron in the hCM frame. 
We also suggest measurement of the asymmetries of the transverse energy 
flow which are simpler and may be calculated reliably. Our predictions for the 
asymmetries of transverse energy flow are the most important contribution of 
this paper. 

1 Large logarithmic corrections and resummation 

The impact parameter resummation formalism that we are applying here de- 
scribes production of nearly massless hadrons in the current fragmentation 
region, where the production rate is the highest. In this region, transverse mo- 
mentum distributions are affected by large logarithmic QCD corrections due to 
radiation of soft and collinear partons. The leading logarithmic contributions 
can be summed through all orders of perturbative QCD [10-12] by applying a 
method originally proposed in [9] for jet production in e + e~ annihilation and 
the Drell-Yan process. 

The spin-averaged cross section for SIDIS in a parity-conserving channel, e.g., 
7* exchange, can be decomposed into a sum of independent contributions from 
four basis functions A p (ip,ip) of the leptonic angular parameters ip,if [13]: 



Here ip is the angle of a hyperbolic rotation (a boost) in Minkowski space; 
it is related to the conventional DIS variable y,byy = P- q/P-£ = 2/(l + 
cosh ■?/>). The angular basis functions are A± — 1 + cosh 2 ^, A2 = —2, A3 = 
— cosy? sinh2-?/>, A 4 = cos2y?sinh 2 '?/>. Of the four structure functions P V, only 
l V and 2 V contribute to the denominator of (1), i.e., the ^-integrated cross 
section. Of these two terms, l V is more singular and it dominates the rate. 
To explore the singular contributions in l V , we introduce a scale qr related 
to the polar angle (9) of the direction of the final hadron in the hCM frame. 
A convenient definition is 

q T = Q^l/x-l exp(-77), (2) 

where 77 is the pseudorapidity of the charged hadron in the hCM frame (de- 
fined with respect to the direction of the momentum g M of 7*). In the limit 
qr — > 0, the structure function l V is dominated by large logarithmic terms; it 

has the form q? 2 Hk=Ms/^) k Em=o v(km) ^(Qt/Q 2 ), where v(km) are some 
generalized functions. To obtain a stable theoretical prediction, these large 
terms must be resummed through all orders of perturbative QCD. The other 
structure functions 2 ' 3 < 4 V are finite at this order; we approximate them by 
fixed-order 0(0:5) expressions. 

In Eq. (1), the numerator of (cosy?) or (cos2y?) depends only on the struc- 
ture function 3 V or A V, respectively. The measurement of (cosy?) or (cos2y?) 
must be combined with good knowledge of the ^-integrated cross section, i.e., 
the denominator of (1), to provide experimental information on the structure 
function 3 V or A V . Thus it is crucial to check whether the theory can re- 
produce the ^-integrated cross section as a function of pr before comparing 
the prediction for (1) to the data. But, on the contrary, as shown in [12], the 
0(a s ) fixed-order cross section is significantly lower than the data from [14] 
in the range of pt relevant to the ZEUS measurements. This difference signals 
the importance of higher-order corrections and undermines the validity of the 
0(a s ) result as a reliable approximation for the numerator of (1). 

On the other hand, the resummation calculation [12] with a proper choice 



of the nonperturbative function yields a much better agreement with the ex- 
perimental data for the (/'-integrated pr-distribution from [14]. One might 
try to improve the theoretical description of the ZEUS data using resumma- 
tion for the denominator of (1). However, the resummation calculation for 
da I (dxdzdQ 2 dp T dip) in the phase space region relevant to the ZEUS data is 
currently not possible, largely because of the uncertainty in the parameteriza- 
tion of the nonperturbative contributions in this region. The resummed struc- 
ture function X V includes a nonperturbative Sudakov factor, which contains 
the effects of the intrinsic transverse momentum of the initial-state parton and 
the nonperturbative fragmentation contributions to the transverse momentum 
of the final-state hadron. Without first determining this nonperturbative fac- 
tor, e.g., from other measurements, it is not possible to make a trustworthy 
theoretical prediction for the denominator of (1) and, hence, these azimuthal 
asymmetries. 

The azimuthal asymmetries measured by ZEUS may also be sensitive to un- 
certainties in the fragmentation to /i ± in the final state. Indeed, the cross 
section in (1) includes convolutions of hard scattering cross sections with frag- 
mentation functions (FFs), integrated over the range 0.2 < z < 1. Although 
the knowledge of FFs is steadily improving [16], there is still some uncertainty 
about their z-dependence and flavor structure for the range of Q relevant to the 
ZEUS measurement. Therefore the most reliable tests of the theory would use 
observables that are not sensitive to the final-state fragmentation. The asym- 
metries {cos rup) would be insensitive to FFs if the dependence on the partonic 
variable z were similar in the hard parts of the numerator and denominator 
of (1), so that the dependence on the FFs would approximately cancel. (We 
denote the parton-level quantities by "~".) It is shown in Appendix B of [11] 
that the partonic structure function l V , which dominates the denominator of 
(1), contains terms proportional to 1/z 2 that increase rapidly as z decreases. 
However, the most singular terms in the partonic structure functions 3 ' 4 V are 
proportional to 1/z. Therefore, the dependence on the FFs does not cancel in 
the azimuthal asymmetries. 

A curious fact appears to support the suggestion that the theoretical predic- 
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Fig. 1. Comparison of the 0(a s ) prediction for the ratio (cos ip) / (cos 2p) with the 
ratio of experimentally measured values of (cosp) and (cos 2p) from [1]. The error 
bars are calculated by adding the statistical errors of (cos p) and (cos 2p) in quadra- 
ture. Systematic errors are not included. The theoretical curve is calculated for 
(x) = 0.022, (Q 2 ) = 750 GeV 2 , using the CTEQ5M1 parton distribution functions 
[15] and fragmentation functions by S. Kretzer from [16]. 

tions for {cosnip) depend significantly on the fragmentation functions. While 
each of the measured asymmetries, (cos ip) and (cos2</>), deviates from the 
0(a s ) prediction, the data actually agree well with the 0(a s ) prediction for 
the ratio (cos ip) /(cos 2ip) , as shown in Fig. 1. The error bars are the statisti- 
cal errors on (cos ip) and (cos 2ip) combined in quadrature; this, however, may 
overestimate the statistical uncertainty if the two errors are correlated. Since 
this ratio depends only on the numerators in Eq. (1), which are less singular 
with respect to z than the denominator, the dependence on the fragmentation 
functions may be nearly canceled in the ratio. The good agreement between 
the 0(a s ) prediction and the experimental data for this ratio supports our 
conjecture that the fragmentation dynamics has a significant impact on the 
individual asymmetries defined in (1). 

Our final remark about the azimuthal asymmetries in (1) is that the Pt (or 
p c ) distributions are not the best observables to separate the perturbative and 
nonperturbative effects. The region where multiple parton radiation effects are 



important is specified by the condition q\IQ 2 <^ 1- But the pt distributions 
are smeared with respect to the qr distributions by an additional factor of 
z, because pt = z qr- Thus the whole observable range of pr is sensitive 
to the resummation effects in the region of qr of the order of several GeV. 
A better way to compare the data to the perturbative QCD prediction is 
to express the azimuthal asymmetries as a function of qr, not pr- Then the 
comparison should be made in the region where the multiple parton radiation 
is unimportant, i.e., for qr/Q ^ 1- 



2 Asymmetry of energy flow 



Next, we describe an alternative test of perturbative QCD, which will fur- 
ther reduce the above theoretical uncertainties: measurement of the azimuthal 
asymmetries of the transverse energy flow. In the hCM frame, the transverse 
energy flow can be written as [17,13,10-12] 

' IL ' ! J2 P V ET (x,Q 2 ,q')A p ^,cp). (3) 



dxdQ 2 dq^dip 



p=i 



Unlike the charged particle multiplicity, the energy flow does not depend on 
the final-state fragmentation. It has been demonstrated [11,12] that a resum- 
mation calculation can provide a good description for the experimental data 
on the (^-integrated -Sp-flow. A new class of azimuthal asymmetries may be 
defined as 

/ rf$ J 2lT cos n<p . J£h-j du? 
{E T cosnip)(q T ) = 2n dEr -f - . (4) 

i d$ io dxdQ^ T d v dl P 

The structure functions p Ve t for the Ex-Row can be derived from the struc- 
ture functions P V for the SIDIS cross section [12]. Similar to the case of the 
particle multiplicities, the asymmetries (Et cos ip) and (Et cos 2ip) receive con- 
tributions from 3 Ve t and 4 Ve t , respectively. But, unlike the previous case, 
the denominator in (4) is approximated well by the resummed E T -Qow. Thus 
these asymmetries can be calculated with greater confidence. 
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Fig. 2. Energy flow asymmetries (Ex cos f)(qx) and (Ex cos 2ip) (qx) for (a) 
x = 0.0047, Q 2 = 33.2 GeV 2 and (b) x = 0.026, Q 2 = 617 GeV 2 . The Figure 
shows predictions from the resummed (solid) and the 0(a s ) (dashed) calculations. 

Figure 2 shows our prediction for the azimuthal asymmetries (Ex cos <p) and 
(£ T cos2<^) as functions of q T for (a) x = 0.0047, Q 2 = 33.2GeV 2 in the left 
plots and (b) x = 0.026, Q 2 = 617GeV 2 in the right plots. The asymmetries 
are shown in g^-bins that are obtained from the experimental pseudorapidity 
bins for the (^-integrated Et-&ow data from Ref. [18]. The upper rc-axis shows 
values of the hCM pseudorapidity r\ that correspond to the values of qr on the 
lower x-axis. For each of the distributions in Fig. 2, the structure functions 
3 Ve t and a Ve t were calculated at leading order in QCD, i.e., 0(a s ). The 
solid and dashed curves, which correspond to the resummed and 0(a s ) results 
respectively, differ because the structure function 1 Ve t in the denominator of 
(4) differs for the two calculations. The resummed (^-integrated .Sr-fiow is 
closer to the data than the fixed-order result, so that the predictions made 
by perturbative QCD for the subleading structure functions 3 Ve t and 4 Ve t 
will be confirmed if the experimental azimuthal asymmetries agree with the 
resummed distributions. 



A recent study [12] shows that in the region qr ~ Q the resummed (^-integrated 
E T -Qow is larger than the 0(a s ) prediction. This explains why the asymme- 
tries for qx ~ Q are smaller for the resummed denominator than for the 0(a s ) 
denominator. In the region qr/Q *C 1, the asymmetries are determined by the 
asymptotic behavior of the fixed-order and resummed partonic structure func- 
tions p Ve t . As qx — > 0, the 0(a s ) structure functions ( vE T )o{a a ), 3 Ve t , and 
a Ve t behave as l/q T , 1/qr and 1, respectively. Thus, asymptotically, the ratios 
3 've t /( vE T )o{a s ) go to zero, although the qr distribution for the asymme- 
try (E T cos<p) is quite large and negative for small, but non- vanishing qT (cf. 
Fig. 2). Resummation of 1 Ve t changes the gr-dependence of the denominator, 
which becomes nonsingular in the limit qx — * 0. Consequently, the asymme- 
try (E T cos<p) with the resummed denominator asymptotically grows as l/qr 
(i.e., in accordance with the asymptotic behavior of 3 Ve t ). Hence neither the 
fixed-order nor the resummed calculation for {Et cos ip) is reliable in the low-gr 
region, so that higher-order or additional nonperturbative contributions must 
be important at qr — > 0. The asymptotic limit for the resummed (E T cos2ip) 
remains finite, with the magnitude shown in Fig. 2. Since the magnitude of 
(Et cos 2^) is predicted not to exceed a few percent, an experimental observa- 
tion of a large asymmetry (E T cos 2ip) at small qT would signal the presence of 
some new hadronic dynamics, e.g., contributions from T-odd structure func- 
tions discussed in [6]. 

Figure 2 shows that the predicted asymmetry (Et cos if) (qx) at qr/Q = 1 is 
about 1-2% for the resummed denominator, while it is about 2-4% for the 
0(a s ) denominator. The asymmetry (E T cos 2ip) (q T ) at qT/Q = 1 is about 
1.5-2% or 3.5-5%, respectively. Both asymmetries are positive for qT ~ Q. 
According to Fig. 2a, the size of the experimental qT bins (converted from the 
77 bins in [18]) for low or intermediate values of Q 2 is small enough to reveal 
the Iow-^t behavior of 3,a Ve t with acceptable accuracy. However, for the high- 
Q 2 events in Fig. 2b, the experimental resolution in qT may be insufficient for 
detailed studies in the low-gr region. Nonetheless, it will still be interesting to 
compare the experimental data to the predictions of perturbative QCD in the 
region qT/Q ~ 1, and to learn about the angular asymmetries at large values 



of Q 2 and x. 

To conclude, we suggest that the azimuthal asymmetry of the energy flow 
should be measured as a function of the scale qr- These measurements would 
test the predictions of the perturbative QCD theory more reliably than the 
measurements of the asymmetries of the charged particle multiplicity. 
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